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ABSTRACT 

We present the first results from a new generation of simulated large sky coverage (^100 square 
degrees) Sunyaev-Zeldovich effect (SZE) cluster surveys using the cosmological adaptive mesh re- 
finement N-body/hydro code Enzo. We have simulated a very large (512 3 h~ 3 Mpc 3 ) volume with 
unprecedented dynamic range. We have generated simulated light cones to match the resolution and 
sensitivity of current and future SZE instruments. Unlike many previous studies of this type, our 
simulation includes unbound gas, where an appreciable fraction of the baryons in the universe reside. 

We have found that cluster line-of-sight overlap may be a significant issue in upcoming single-dish 
SZE surveys. Smaller beam surveys (~1') have more than one massive cluster within a beam diameter 
5-10% of the time, and a larger beam experiment like Planck has multiple clusters per beam 60% of 
the time. We explore the contribution of unresolved halos and unbound gas to the SZE signature at 
the maximum decrement. We find that there is a contribution from gas outside clusters of ^16% per 
object on average for upcoming surveys. This adds both bias and scatter to the deduced value of the 
integrated SZE, adding difficulty in accurately calibrating a cluster Y-M relationship. 

Finally, we find that in images where objects with M > 5 x 10 13 M have had their SZE signatures 
removed, roughly a third of the total SZE flux still remains. This gas exists at least partially in 
the Warm Hot Intergalactic Medium (WHIM), and will possibly be detectable with the upcoming 
generation of SZE surveys. 

Subject headings: cosmology: theory-galaxies:clusters:general-cosmology:observations- 
hydrodynamics-methods:numerical-cosmology:cosmic microwave background 



1. INTRODUCTION 

Clusters of galaxies form from the highest peaks in the 
primordial spectrum of density perturbations generated 
by inflation in the early universe. They are the most mas- 
sive virialized structures in the universe, and as such are 
rare objects. The number density of galaxy clusters as a 
function of mass and redshift is strongly dependent on a 
number of cosmological parameters. In particular, count- 
ing the abundance of clusters above some lower mass 
limit as a function of cluster redshift places constraints 
on Of,, Qm-i cs , and the dark energy equation of state 
param eter, w (W ang fe SteinhardtJ Il998t lHaiman et al.l 

Observational measurement of the cluster abundance 
over large sky areas and redshift range is required to 
generate cosmological parameter constraints which are 
complementary to constraints from the cosmic microwave 
background (CMB), type la supernovae, Big Bang nucle- 
osynthesis (B BN), the Ly man-q forest, and galaxy red- 
shift surveys (lTozzill2006li . Cluster survey yields depend 
on the value of the minimum cluster luminosity probed 
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as a function of redshift, the scaling between cluster lu- 
minosity and mass, the growth function of structure, and 
the redshift evolution of the comoving volume element, 
all of which depend on a complex com bination of cosmo- 
logical parameters dRosati et al.ll2002f) and intracluster 
medium (ICM) physics (|Evrardll2004D ' 

1.1. Sunyaev-Zel'dovich Effect Surveys 

Th e Sunyaev-Zel'dovich Effect (jSunvaev fc Zeldovichl 
Il972f ) is a process by which the hot electrons trapped 
in the large dark matter potential wells of clusters in- 
verse Compton scatter CMB photons to higher energy, 
resulting in a low frequency (< 218 GHz) decrement, and 
a corresponding high frequency (> 218 GHz) increment 
in the CMB on the angular scale of the cluster. The 
strength of the SZE decrement/increment is character- 
ized by the Compton y parameter, which results from 
the line of sight integral of the thermal pressure 

y = \ a T n e ^—dl, (1) 
J m e c z 

where n e is the electron density and T is the gas electron 
temperature. We also define the integrated Compton y 
parameter as 

Y = JydA (2) 

which is the integration over the area subtended by a 
circle corresponding to some relevant cluster physical ra- 
dius. The observed temperature fluctuation correspond- 
ing to a give value of y in a given frequency band is 
AT 



2 



where 

g(x) = (*£^k - 4) [1 + S„ e (x,T e )], (4) 

x = hv/kTc mb and Ssze is a rela tivistic correction as 
described in lltoh fc Nozawal (|2004h . For our purposes in 
this paper we have neglected the relativistic correction, 
but will explore it in future work. This correction is small 
(less than 1% at the maximum decrement frequency) for 
clusters with T < 10 keV. 

The SZE is particularly useful in cos mological stud- 
ies due to its nea r redshift independence (Rcphacli] |1995l : 
lBirkinshawlll999t ICarlstrom et al.ll2002ft . Therefore, ob- 
servations of clusters are not limited to low redshift as in 
the X-ray, but can extend to as high as z ~ 2, where the 
number of massive clusters becomes small. An additional 
consequence of the redshift independence of SZE sur- 
veys is that a flux-li mited survey is also approximately a 
mass limited survey (jRephacli 19951 lHaiman et al.|[200ll ). 
These two properties make SZE surveys uniquely valu- 
able for cluster abundance counts and determination of 
the cluster mass function with redshift, provided one can 
obtain independent optical redshifts for the galaxies in 
the identified objects. The near redshift independence of 
the SZE creates unique complications for large area sur- 
veys which do not seriously affect other types of surveys 
(e.g., optical and X-ray). In particular, the contribution 
to the sky signal from both low mass and distant halos, 
as well as unbound gas, may be a significant source of 
confusion. 

There are several upcoming millimeter wavelength 
cluster surveys with new telescopes including the Ata- 
cama Pathfind er Experiment Su nyaev-Zeldovich survey 
(APEX-SZ) (iGiisten et all 12 006). the South Pole Tele- 
scope (SPT) (jRuhl et al.ll2004D. and th e Atacama Cos- 
mology Telescope (ACT) (jfowlerl I2004T) . in add ition to 
the space-based Planck Surveyor (|Murdinll2000| ). which 
will conduct large blind surveys of clusters using the SZE. 
Table [1] shows the values of the relevant instrumental 
characteristics for the single-dish survey telescopes we 
have taken from the literature and used in the follow- 
ing analysis. For this study, we limit ourselves to re- 
sults using a single band, ^144 GHz, where the SZE 
decrement is maximal, though these results will general- 
ize to other bands and to multiwavelength studies, since 
the SZE from all gas will have roughly the same spec- 
tral signature (modified slightly by relativistic effects). 
However, for removing contaminating signals (e.g. radio 
point sources, CMB) multiwavelength coverage will be 
very desirable. 

The variation in survey characteristics for these instru- 
ments has important consequences for cluster surveys. 
For example, the distribution of sources detected as a 
function of cluster redshift should be different for each 
survey. This is because although the SZE surface bright- 
ness does not diminish with distance, the angular size of 
the objects does vary with redshift, and may be larger or 
smaller than the instrument beam for any given cluster. 
This selection effect is modified also by the volume sam- 
pled as a function of redshift in a fixed angular field and 
the growth rate of structure in an ACDM universe. It is 
important to understand these selection effects in order 
to constrain cosmology. We must be able to determine 
the correct distribution of clusters as a function of red- 



TABLE 1 

Characteristics of Upcoming SZE Surveys 



Survey 


Angular Coverage 


Beam Size (~144 GHz) 


RMS per beam 


APEX-SZ 


TBD 


1.0' 


10 AiK 


SPT 


4000 deg 2 


1.0' 


10 //K 


ACT 


100 deg 2 


1.7' 


2 nK 


Planck 


All-Sky 


7.1' 


6.0 fiK 



shift from the surveys or systematic errors in estimated 
cosmological parameters will result. 

There are also several centimeter wave («30 GHz) 
interferometers that are performing surveys of the 
SZE, such as the Arcminute Microkelvin Imager (AMI) 
dHollerll2000h an d the Sunyaev-Zel'dovich Array (SZA) 
(|Loh et al.l l2005f ). SZE surveys have the potential to 
strongly constrain the w parameter for dark energy, since 
they sample clusters to large redshift. 

Determining the abundance and distribution in mass 
and redshift of massive galaxy clusters from observables 
is a challenging exercise. In the realm of cluster abun- 
dance counts, one needs to know with high precision the 
mass range of clusters probed in a flux-limited survey as a 
function of cluster redshift. This determination depends 
on detailed knowledge of the scaling between mass and 
light in clusters. It is critical to understand how cluster 
observables correlate with cluster total mass in order to 
use clusters of galaxies as precision cosmological tools. 

1.2. The Role of Simulations in Understanding Cluster 

Surveys 

Recent res ults indicate that high resolution N-body 
simulations (I Warren et aH 120061 : iHeitmann et alj 120061 : 
I Reed et all l2005f) generate mass functions which differ 
significantly from the Press-Schechte r result and also 
from the s ubsequent modificati ons of ISheth fc TormenI 
(|1999f) and Uenkins et all (|2001fh Since there is strong 
evidence that purely analytic methods are insufficient, 
"precision cosmology" requires the use of numerical sim- 
ulations. In other words, in order to make predictions 
which match the observed cluster population to percent- 
level precision, analytic methods are inadequate. 

The output of numerical simulations of clusters can be 
compared to the observed cluster mass function. This 
comparison is non-trivial, however, due to the uncer- 
tain nature of the conversion between observable quanti- 
ties (e.g., X-ray luminosity, SZE Compton y parameter, 
lensing shear) and cluster total mass. Observations of 
the cluster gas typically depend on the detailed proper- 
ties of the hot baryons in clusters. It has been shown 
by our group and others that cluster observables have a 
strong dependence on the baryoni c physics in the ICM 
(jHallman et al.l 120061 : iNagail l2006h . and are subject to 
an array of uncertainties. However, the SZE signal in- 
tegrated over the projected cluster area (as defined in 
Equation 2) inside rsoo is unique in that the general scal- 
ing with ma ss is relatively in d ependent of the assumed 
gas physics (jMotl et all 120051 : lNagaill2006l ). While the 
normalization of the Y-M relation has some dependence 
on ICM physics, th e slope and tightness of the correla- 
tion are unaffected (Nagai 2006). Most recent simulated 
light cone calculations use the dark matter mass function 
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generated by large N-body simulations like the Hubble 
Volume simulation or the Millennium run, with "painted 
on" baryons to generate moc k surveys (e.g-. lEvrard et al] 
I2002t iGeisbusch et al.ll2005[) . These studies typically as- 
sume the gas is isothermal and in hydrostatic equilibrium 
with the dark matter potential. 

Both simu l ations including r elevant physics (e.g. , 
iWhitd 12001 iRasia et alJ [200l lHallman et all l2006h 
and high resolution X-ray observations of galaxy clus- 
ters (e.g.. IVikhlinin et"al]|2005t iMarkevitch et all 120021 : 
lHallman fe Markevitchl l2004f ) suggest that many clus- 
ters depart strongly from both equilibrium and isother- 
mality. These deviations can have a strong impact on 
both the observable and derived properties of clusters. 
Thus, in order to properly simulate sky surveys, it is 
critically important to self-consistently include baryons 
in numerical simulatio ns. While some work has been 
done in this area (e.g., jSpringel et al.ll200ll : I White et al.1 
120021 iRoncarelli et al.l l2006l 120071 ). the largest volumes 
simulated were small (~100-200 h _1 Mpc), and only suf- 
ficient to generate synthetic light cones of roughly 1-4 
deg 2 . The simulation performed for this study models a 
significantly larger physical volume than previous efforts, 
has a higher peak physical resolution, and fully incorpo- 
rates baryons self-consistently. This allows us to perform 
much larger synthetic surveys (^100 deg 2 ) than could be 
done with earlier N-body/hydro simulations. 

Cosmological N-body/hydro simulations have ad- 
vanced significantly in the last decade, such that the 
simulation output now compares quite well to observa- 
tions of galaxy clusters (|Springel et al.ll2006h . In partic- 
ular, our group and others have shown that there is good 
agreement in simulated and observed scaling relations be- 
tween bulk cluster ICM properties (e.g., clust er mass, X- 
ray lu minosity, X-ray spectral temperature) (jMotl et al.l 
2004). There remain important differences, particularly 
for lower mass clusters, which indicates the need for a 
better understanding of the details of baryonic cluster 
physics. The advance of realism in simulations is a result 
of diligent, iterative efforts by various groups of investi- 
gators to directly compare simulations to observations. 
This study uses a large volume high resolution adiabatic 
simulation, and serves as a template for more complex 
runs involving additional non-gravitational physical pro- 
cesses which are likely important to accurately modeling 
cluster surveys. These results should be relatively ro- 
bust in any case, as it has been shown that SZE survey 
yields are relatively independent of cluster physics details 
flWhite et al.ll2002T) . 

1.3. Modeling SZE Surveys 

A variety of approaches have been taken to model SZE 
surveys. Most involve either semi-analytic prescriptions 
or N-body simulations where the gas is adde d in a post 
processing step (e.g., ISchulz fe White! |2003). As dis- 
cussed in the previous section, there are limitations to 
these methods, particularly in the assumptions of hydro- 
static equilibrium and isothermality. Some studies (e.g., 
IWhite et al]l2002HHolder et al.ll2007ft have discussed the 
contribution of unresolved clusters/groups to the signa- 
ture of detected clusters. The presence of gas outside the 
cluster virial radii in low density unbound structures such 
as filaments is potentially also quite important. This gas 
is completely absent in non-hydrodynamic treatments of 



the simulation volume, appears naturally in our calcula- 
tion, and is exp ected to contain 40-50 % of the baryons 
in the universe (|Cen fc Ostrikerll2006h . Since the SZE 
does not diminish with distance, and results from a line 
of sight integral of the gas pressure, the sum of all the 
flux from unbound gas could be a significant contributor 
to the total flux in any cluster detection. 

In this paper, we have examined the contribution of 
line-of-sight baryonic gas to the expected SZE signal of 
simulated clusters. We have stacked a (512 h _1 Mpc) 3 
volume (comoving) adaptive mesh refinement (AMR) N- 
body/hydro simulation to generate a survey of a light 
cone subtending 100 square degrees on the sky. Our 
simulated survey covers a larger sky area, and at higher 
angular resolution, than any previous N-body + hydro 
simulated survey. 

There is extensive work in the literature o n cluster de- 
tection algorithms for SZE surveys, (e.g.. iDiego et al.l 
|2002t iHerranz et all 12003 iHobsonfe McLachlanl 120031 : 
ISchafer et ail l2006t IMelin et all I2006D . These methods 
involve various techniques designed to spatially filter out 
the primary CMB anisotropies (so called matched filter- 
ing), wavelet techniq ues, and application of p ublic tools 
such as SExtractor 6 (|Bertin fc Arnoutsl IT996) . We find 
the existing work to be quite detailed, and do not in- 
troduce new algorithms of this type here. Indeed, it is 
important to step back from the analyses which have 
attempted to include all the relevant contaminants and 
instrumental effects and explore the intrinsic difficulties 
resulting from the cluster population as projected on the 
sky. There is a limiting precision one can expect from 
cluster surveys irrespective of the ability to remove in- 
strumental effects, point source confusion and sky back- 
grounds. This limit results from the possibly irreducible 
confusion due to clusters, groups, lower mass halos and 
unbound gas, all of which contribute to the SZE signal 
with a nearly identical spectral signature. This study 
examines these effects with a more realistic cosmological 
calculation than has typically been done, including the 
full complement of baryons expected in the real universe. 

We explore the intrinsic limitations of SZE surveys by 
comparing and characterizing the contribution of unre- 
solved halos and unbound filamentary gas to the cluster 
signal in samples that might result from upcoming sur- 
veys using a full hydro/N-body simulation of the cluster 
sky. In future work, we will model backgrounds and in- 
strumental charact eristics as has been done in the litera- 
ture recently (e.g iSehgal et al.l l2006t IMelin et al.l 120061 : 
ISchafer et al.ll2006h with a focus on techniques for accu- 
rately extracting cluster properties and abundance. 

We discuss our methodology of simulating these large 
surveys in Section 2, analytic predictions of SZE observ- 
ables in Section 3, present results in Section 4, discussion 
in Section 5, and summarize our work in Section 6. 

2. METHODOLOGY 

2.1. The Enzo Code 

'Enzo' 7 is a publicly available, extensively tested 
adaptive mesh refinement co smology code developed b y 
Greg Bryan and colleagues dBrvan fc Normanl fl997al lbl: 
INorman fc B7va^[l999l ' lO'Shea et al.ll2004L 120051) . The 

6 http:/ /tcrapix.iap.fr/rubriquc.php?id_rubrique=91/ 

7 http:/ /lea. ucsd.edu/portal/software/enzo/ 
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specifics of the Enzo code are described in detail in these 
papers (and references therein), but we present a brief 
description here for clarity. 

The E nzo code couples a n N-b o dy particle-mesh (PM) 
solve r (jEfstathiou et al.l 119851 : iHocknev fc Eastwood! 
1 1988T ) used to follow the evolution of a collisionless dark 
matter component with an Eulerian AMR met hod for 
ideal gas dynamics by iBerger fc Colellal (|1989| ). which 
allows high dynamic range in gravitational physics and 
hydrodynamics in an expanding universe. This AMR 
method (referred to as structured AMR) utilizes an adap- 
tive hierarchy of grid patches at varying levels of res- 
olution. Each rectangular grid patch (referred to as a 
"grid") covers some region of space in its parent grid 
which requires higher resolution, and can itself become 
the parent grid to an even more highly resolved child 
grid. Enzo's implementation of structured AMR places 
no fundamental restrictions on the number of grids at 
a given level of refinement, or on the number of levels 
of refinement. However, owing to limited computational 
resources it is practical to institute a maximum level of 
refinement, £ m ax- Additionally, the Enzo AMR imple- 
mentation allows arbitrary integer ratios of parent and 
child grid resolution, though in general for cosmological 
simulations (including the work described in this paper) 
a refinement ratio of 2 is used. 

Since the addition of more highly refined grids is adap- 
tive, the conditions for refinement must be specified. 
In Enzo, the criteria for refinement can be set by the 
user to be a combination of any or all of the following: 
baryon or dark matter overdensity threshold, minimum 
resolution of the local Jeans length, local density gra- 
dients, local pressure gradients, local energy gradients, 
shocks, and cooling time. A cell reaching any or all of 
the user-specified criteria will then be flagged for refine- 
ment. Once all cells of a given level have been flagged, 
rectangular solid boundaries are determined which min- 
imally encompass them. A refined grid patch is then 
introduced within each such bounding volume, and the 
results are interpolated to a higher level of resolution. 

In Enzo, resolution of the equations being solved is 
adaptive in time as well as in space. The timestep in 
Enzo is satisfied on a level-by-level basis by finding the 
largest timestep such that the Courant condition (and 
an analogous condition for the dark matter particles) is 
satisfied by every cell on that level. All cells on a given 
level are advanced using the same timestep. Once a level 
L has been advanced in time At l , all grids at level L + 1 
are advanced, using the same criteria for timestep calcu- 
lations described above, until they reach the same phys- 
ical time as the grids at level L. At this point grids at 
level L + 1 exchange baryon flux information with their 
parent grids, providing a more accurate solution on level 
L. Cells at level L + 1 are then examined to see if they 
should be refined or de-refined, and the entire grid hier- 
archy is rebuilt at that level (including all more highly 
refined levels). The timestepping and hierarchy rebuild- 
ing process is repeated recursively on every level to the 
maximum existing grid level in the simulation. 

Two different hydrodynamic methods are implemented 
in Enzo: the piecewise parabolic method (PPM) 
(I Woodward fc Colellalfl98l . w hich was extended to cos- 
mology by iBrvan et al.l (|1995f ). and the hydrodynamic 
method used in the ZEUS magnctohydrodynamics code 



(| Stone fc Normanl ll992alTbT ) . We direct the interested 
reader to the papers describing both of these methods 
for more information, and note that PPM is the preferred 
choice of hydro method since it is higher-order-accurate 
and is based on a technique that does not require artifi- 
cial viscosity, which smoothes shocks and can smear out 
features in the hydrodynamic flow. 

2.2. Simulation Setup and Analysis 

The simulation discussed in this paper is set up as fol- 
lows. We initialize our calculation at z — 99 assuming a 
cosmological model with f2 m = 0.3, f2& = 0.04, Hcdm — 
0.26, n A = 0.7, h = 0. 7 (in units of 100 km/s /Mpc), 
(is = 0.9, and using an lEisenstein fc Hul (|1999t ) power 
spectrum with a spectral index of n — 1. The simulation 
is of a volume of the universe 512 h _1 Mpc (comoving) 
on a side with a 512 3 root grid. The dark matter parti- 
cle mass is 7.228 x 10 10 h _1 M Q and the mean baryon 
mass resolution is 1.112 x 10 10 h _1 M©. The simulation 
was then evolved to z = with a maximum of 7 levels of 
adaptive mesh refinement (a maximum spatial resolution 
of 7.8 h^ 1 comoving kpc), refining on dark matter and 
baryon overdensities of 8.0 (to ensure an approximately 
Lagrangian mass resolution in baryonic structures) . The 
equations of hydrodynamics were solved with the Piece- 
wise Parabolic Method (PPM) using the dual energy for- 
malism. The entire grid hierarchy (including both parti- 
cle and baryon information) was written out at regular 
intervals, and in particular, data was output at intervals 
of Az = 0.25 between z = 3 and z = 2.5 (inclusive), and 
Az = 0.1 between z = 2.5 and z = 0.1 (inclusive). 

Analysis was performed on every data output between 
z = 3 and z = 0.1 in an identical way. Th e HOP halo- 
finding algorithm (jEisenstein fc Hutl Il998l ) was applied 
to the dark matter particle distribution to produce a 
dark matter halo catalog. Spherically-averaged, mass- 
weighted radial profiles of various baryonic and dark mat- 
ter quantities including density, temperature, and pres- 
sure were then generated for every halo in the catalog 
with an estimated halo mass greater than 4 x 10 13 M . 
These radial profiles were used to calculate more accu- 
rate virial masses and radii as well as an estimate for the 
Compton y parameter as a function of impact parameter 
on the halo. Projections of the integrated Compton y 
value along the line of sight were created for each of the 
three axes along the simulation volume, with two projec- 
tions per axis - one of the front half of the simulation 
volume, and one along the back half. Each projection 
has an approximate depth of Az = 0.1. These projec- 
tions have a resolution of 2048 pixels on a side. 

2.3. Generation of the "Santa Fe" Light Cone 

Mock SZE observations of a 10° x 10° patch of sky 
are generated by stacking projections from the simula- 
tion discussed in Section 12.21 These "light cones" are 
created by stacking projections of the Compton y pa- 
rameter at each redshift output. At each redshift, the 
projection is chosen to be along a random axis, and has 
been randomly shifted in space such that the positions 
of large scale structure is uncorrelated. Additionally, the 
projections have been rescaled to the resolution of the 
light cone, which is 2048 pixels per side, or a resolution 
of 17.58 arc seconds per pixel. This scaling may involve 



5 



tiling (for redshifts where 512 h _1 Mpc comoving corre- 
sponds to less than 10° on the sky) or interpolating (for 
redshifts where 512 h _1 Mpc comoving corresponds to 
more than 10° on the sky). Secondary maps are created 
which include only the Compton y parameter contributed 
by gas within the virial radius of halos with masses above 
5 x 10 13 Mq, and only gas outside of the virial radius of 
these objects. 200 of these mock "light cones" at this size 
and resolution are created using different random seeds. 
These light cones (named "Santa Fe" light cones due to 
the location where the project was conceived) have angu- 
lar resolution which is significantly higher than any cur- 
rent or proposed SZE observational campaign. The goal 
of this analysis is not to determine an optimal method for 
source finding, but to determine the contamination from 
unresolved halos and unbound gas for a simple method. 

3. ANALYTIC PREDICTIONS FOR SZE OBSERVABLES 

Here we describe some of the theory behind the use 
of SZE cluster observations in constraining cosmology. 
Though we are aware that these types of analytic calcu- 
lations have been performed previously, we show them 
here to motivate not just the current analysis, but that 
which will be performed for subsequent papers in this 
series. 

One of the most useful methods for retrieving cosmo- 
logical information from SZE observations of galaxy clus- 
ters is by the calculation of galaxy cluster counts as a 
function of redshift. The number of galaxy clusters above 
some given minimum mass M m i n (z) in a redshift bin of 
width dz and solid angle dp. can be defined using the 
Press-Schechter formalism (|Press fe S chcchterl ll974f ) as 



dN 
dzdfl 



dV 
dzdfl 



dM^{M,z) (5) 



where dV/ dzd£l is the cosmological comoving volume el- 
ement at a given redshift, and -j^dM is the comoving 
halo number density as a fu nction of mass and re dshift. 
The latter is expressed as by I Jenkins et al.l (|2001l ) as 



dn 
dm 



(M,z) 



-0.315 



Po 1 du M 



M om dm 
exp [-|0.61 - log{D(z)a M )\ 3 - s ] (6) 

where o~m is the RMS density fluctuation, computed on 
mass scale M from th e z = linear power spectrum 
(|Eisenstein fe Hul [T999), p is the mean matter density 
of the universe, defined as po = fl m Pc (with p c being the 
cosmological critical density, defined as p c = 3Hq /8ttG) , 
and D(z) is the linear growth function, given by this fit- 
ting function: 

1 5Q m (*) 
V ; 1+2 2 

l{z fn n A (*) + [i + + ^]p7) 

(|Carroll et al.lfl992t ). with fl m (z) and fl\(z) defined as: 

n m (z) = n mfi (i + zfE- 2 (z) (8) 

and 

n A (z) = n Afl E- 2 (z) (9) 

where Q m ,Q and ^a,o are the density of matter and cos- 
mological constant at the present day, expressed in units 



of the critical density. The cosmological volume element 
is given by: 



dV 
dzdQ 



(*) = 



Ha 



(l + z) 2 D* 
E(z) 



(10) 



where Da(z) is the angular diameter distance as a func- 
tion of redshift, c is the speed of light, Ho is the Hubble 
parameter at z = 0, and E(z) is given by: 

e 2 {z) = n mfi (i + z) 3 + n A (11) 

in a flat universe with a cosmological constant (jPeeblesI 
1993). The RMS amplitude of the density fluctuations 
as function of mass at any redshift, as smoothed by a 
spherically symmetric window function with a charac- 
teristic comoving radius R, can be computed from the 
matter power spectrum using the relation: 



a 2 (M,z) 



f^P(k,z)\W R (k)\ 2 (12) 



where Wn(k) is the Fourier transform of the real-space 
top hat smoothing function: 



k 3 R 3 



i(kR) - kRcos(kR)] (13) 



The radius R is calculated for a given mass by using 
the relation M = ^nR 3 Q m p c , and a(M,z) is normal- 
ized to as, defined as the RMS density fluctuation when 
smoothed by a sphere with a comoving radius of 8 h _1 
Mpc at z = 0, using observations of large-scale structure 
or the cosmic microwave background. The matter power 
spectrum is expressed as: 



2^2 



P(k,z) 



ck 

#0 



3+n 



T 2 (k) 



DM 
D 2 (0) 



(14) 



where T(k) is the matter transfer function which de- 
scribes the way in which the processing of the initial spec- 
trum of matter density fluctuat ions during the radiation- 
dominated era ([Peebles! 119931 ) and D(z) is the fitting 
function for the linear growth function, as given in Equa- 
tion [7] In the calculations discussed in t his paper, we use 
the tr ansfer function T{k) provided by lEisenstein fe Hul 
(|1999h . 

Figure Q] shows the number of galaxy clusters per 
square degree as a function of redshift with M tot > 1 X 
10 14 h" 1 M Q in the WMAP Year III "most favored" cos- 
mology (Q m = 0.268, erg = 0.776) and in the cosmology 
used in the simulation in this paper (fi m = 0.3, as = 0.9). 
Due to the higher il m and as, significantly more galaxy 
clusters are expected to be seen in this simulation than 
one would expect given the WMAP Year III result. 

Figure shows the number of galaxy clusters per 
square degree as a function of redshift with M to t > 
1 x 10 14 h _1 M Q for a variety of cosmological models. 
Panel (a) shows a sequence of cosmologies where all pa- 
rameters except as are held constant, using the same 
cosmological parameters as in the simulation described 
in this paper (fi TO = 0.3, Cl b = 0.04, £l A = 0.7, h = 0.7, 
n = 1), but varying as from 0.6 to 1.2 in steps of 0.1 (bot- 
tom to top lines). Panel (b) shows a sequence of cosmolo- 
gies where all parameters except f2 m and £l\ are held con- 
stant using the same cosmological parameters as in the 
simulation described in this paper (fib = 0.04, h = 0.7, 




Fig. 1. — Number of galaxy clusters per square degree as a func- 
tion of redshift with Mtot > 1 X 10 14 h" 1 M in the WMAP Year 
III "most favored" cosmology (solid line; Q m = 0.268, <rg = 0.776) 
and the cosmology used in the simulation in this paper (dashed 
line; f2 m = 0.3, erg = 0.9). The distribution peaks at z ~ 0.55 for 
the WMAP Year III cosmology and at z ~ 0.64 for the cosmology 
used in the simulations in this paper. 

n = 1, o"8 = 0.9), but varying tt m from 0.2 to 0.4 in steps 
of 0.05 (bottom to top lines), and keeping J7a = 1 — Q m - 
It is interesting to note that varying cr 8 while holding 
all other cosmological parmaters results in a change in 
both the overall number of halos (a factor of more than 
3 between WMAP and the cosmology in our simulation) 
and the redshift at which the distribution peaks . Both 
effects can be explained by examining the term in the 
exponent in equation [6] The comoving number density 
of halos is maximized when 0.61 — log(D(z)au) = 0, 
or D{z) x om — 4.07. For a given set of cosmological 
parameters, increasing as increases om for a given mass 
value, and thus maximizes the number density of halos at 
a given mass at a smaller linear growth factor (or, more 
intuitively, a higher redshift). Variation in tt rn and fl\ 
while holding all other parameters constant effectively 
results in a change in the normalization of the overall 
halo number density, while keeping the redshift at which 
this distribution peaks roughly constant. 

These results are of great interest for upcoming SZE 
surveys, which will sample clusters to relatively low mass 
at high redshifts (z = 0.5-1.0) compared to optical and 
X-ray surveys. This redshift range is where we should 
expect the largest difference between a low a$ cosmology 
preferred by WMAP and a higher value typically used in 
cosmological simulations. The large difference between 
the abundance of clusters in these different cosmologies 
should lead to large differences in the number of identified 
clusters in surveys. This means that very early on in 
any SZE survey it should become fairly obvious which 
cosmology is preferred by cluster observations. This will 
be an even stronger constraint when optical follow up 
observations are used to determine redshifts, which will 
break the degeneracy between VLm and erg. 
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Fig. 2. — Number of galaxy clusters per square degree as a func- 
tion of redshift with Mtot > 1 X 10 14 h -4 Mq for a variety of 
cosmologies. Panel (a): All parameters except <rg are held con- 
stant, using the same cosmological parameters as in the simulation 
described in this paper (Q m = 0.3, O;, = 0.04, Hyv = 0.7, h = 0.7, 
n = 1), but varying as from 0.6 to 1.2 in steps of 0.1 (bottom to 
top lines). Panel (b): All parameters except Q m and Qa are held 
constant using the same cosmological parameters as in the simula- 
tion described in this paper (f2j, = 0.04, h = 0.7, n = 1, <rg = 0.9), 
but varying Q m from 0.2 to 0.4 in steps of 0.05 (bottom to top 
lines) and keeping Q\ = 1 — Q m . 

4. RESULTS 

4.1. Simulated Mass Function 

One of the most basic tests of the correctness of a cos- 
mological simulation is whether or not it can match the 
predicted halo mass function for a given cosmology. This 
is particularly important in the context of creating sim- 
ulated sky maps for cosmological surveys of any kind, 
given that the number of halos and their redshift distri- 
bution is the most basic test of cosmological properties. 
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Fig. 3. — Cumulative number density at several redshifts. 
Solid: z = 0.1. Short-dashed: z = 1.0. Long-dashed: 2 = 2. 
Thick lines are simulation results calculated using the halo virial 
masses (dm+ gas) and the thin lines are the Warren fitting function 
llWarren et al.ll2006h . 

Figure [3] shows the cumulative number density of cosmo- 
logical halos as a function of mass for several redshifts 
which span the range of interest for the topics discussed 
in this work (z = 0.1 — 3). The halos were found as de- 
scribed in Section [221 and the masses used are the total 
halo virial mass (including both baryons and dark mat- 
ter), rather tha n the mass returned by t he Hop halo find- 
ing algorithm (|Eisenstein fc HutJ|1998f ). We also show 
the fitting f unction for cumulati ve halo number density 
obtained bv lWarren et al. (2006) as a reference. At low 
redshifts [z = 0.1 — 1) the mass function of halos from 
the simulation agrees quite well with the fitting function 
over the mass range of interest. This is encouraging, as 
the bulk of galaxy clusters in the universe are at these 
low redshifts (as shown in Figure [T]). At higher redshifts 
(z = 2), the fitting function and halo mass function only 
agree at the highest masses (Mhaio > 10 14 /i -1 M Q ). This 
is to be expected: grid-based codes, including adaptive 
mesh codes, tend to suppress low-mass halo formation, 
particularly at hig h redshift, as has be e n seen in recent 
code comparisons (O' Shea et al.l l2005t iHeitmann et al.l 
2005). This suppression in our study only exists due 
to the choice of simulation setup. In order to model such 
a large physical volume with both N-body + hydrody- 
namics, we must sacrifice mass resolution due to com- 
putational effort concerns. Given that the cosmological 
surveys of interest will only be sensitive to halos in the 
mass range where this simulation agrees well with the 
Warren et al. fit, and the relative paucity of these halos 
at z > 2, there is little cause for concern. 

4.2. SZE Angular Power Spectrum 

We have calculated angular power spectra for the 200 
survey images to determine the cosmic variance in this 
field, and the general form of the power spectrum. The 
calculation involves determining the power spectral den- 
sity of the image as a function of image scale. Addition- 



ally, we have performed the same calculation on maps 
smoothed with the beam size at 144 GHz for each of 
four upcoming SZE survey instruments. We have modi- 
fied the images to model the limitations of these instru- 
ments in a simple way. First, for each instrument, we 
have gaussian smoothed the image with the FWHM of 
the beam at a wavelength (frequency) of 2.1mm (144 
GHz) corresponding to the maximum SZE decrement. 
All the upcoming SZE single-dish surveys will have the 
capability of operating at this wavelength. We have also 
rcbinncd each image such that the beam diameter is rep- 
resented by at least two pixels in order to be Nyquist 
sampled. A "background" in each case is generated by 
adding a gaussian distributed variation with a FWHM 
equal to each survey's stated limiting sensitivity. 

The result of our analysis of the power spectrum of 
the simulated SZE surveys is shown in Figure SJ The 
solid lines are the mean values of the power from the 200 
stacking realizations of the light cone, and the dotted 
lines indicate the range in which 90% of our simulated 
light cone power spectra fall. The color indicates which 
survey's characteristics were used to generate the result. 

The location of the peak of each curve is a function 
of the resolving power of the survey, as is its amplitude, 
in that the power at the smallest angular scale for each 
survey is different. Th e cyan line (la b eled S WH01) is 
the result obtained by ISpringel et al.l (12001 ■ with a 1 
degree angular scale light cone generated from an SPH 
N-body/hydro simulation run with cosmological param- 
eters matching our simulation. The raw power spectrum 
is generated with no smoothing from the raw light cone 
survey images. The SWH01 result is different from ours, 
possibly a result of lower spatial resolution in the sim- 
ulation (resulting in a deficit in the small scale angular 
power) and had a smaller angular field (1 square degree) 
than ours, leading to a higher i cutoff (i?=400) than we 
show (£=70). As ha s been shown p r eviou sly, and specif- 
ically in Figure 4 of ISpringel et al. (2001), analytic pre- 
dictions of the SZE power vary widely depending on the 
calculation method. They neither agree with the sim- 
ulation results, nor in many cases with each other. It 
seems that the analytic result has not converged, and so 
to avoid confusion we do not plot it in this work. 

It is clear that each survey will sample a slightly dif- 
ferent range of scales, though they obviously all are able 
to measure the large angular scale power. Each power 
spectrum peaks and turns over where the angular scale 
of the beam begins to limit the measurement at high 
(,. The power spectrum is not very sensitive to non- 
gravitational gas physics at low multipole numbers, but 
prima rily is sensitive to cosmology, particul arly the value 
of cr 8 dWhite et al.l[200l : iHolder et al.ll200l . The addi- 
tion of non-gravitational physics d oes impact th e smal l 
scale power however, for example Hol der et al.l (|2007l ) 
show that preheating results in reduced small scale power 
in their simulated images. 

What is also of interest in this analysis is the size of 
the variance shown by the 90% range error bars. On 
100 square degree patches of the sky at t < 2000 the 
deduced power can be different by factors of 5-8. This 
indicates that the power spectrum can be quite different 
from one area of the sky to another, and clearly requires 
greater sky coverage to be well constrained. The cosmic 
variance range does not become very small until i > a 
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Fig. 4. — Upper Left: Light cone survey image of 100 square degrees modified to simulate the beam size and limiting sensitivity of the 
Planck Surveyor all-sky survey at 144 GHz (7.1', 6/jK). Upper right: Angular power spectra generated from these images. Lower left: Light 
cone image modified to simulate APEX-SZ/SPT survey characteristics at ~144 GHz (1.0', 10/iK). Lower right: Light cone image modified 
to simulate ACT survey (1.7', 2/xK). Numbers in parentheses indicate (beam size, survey sensitivity /beam) at 144 GHz. 



few thousand. 

4.3. SZE Source Identification 

To identify objects in the light cone images, we simply 
locate the projected clusters from the three-dimensional 
halo finding in the image plane. Since we design the shift- 
ing and stacking strategy, it is trivial to determine the 
image plane location of each cluster in each redshift slice 
of the light cone. We also have calculated the spherically 
averaged radial profiles projected into the image plane 
of the Compton y-parameter. For each cluster with M 
> 10 14 M Q , we can then calculate the integrated value 
of Y (when comparisons to analytic results are not per- 
formed, we use true masses on the simulation grid, with- 
out the h^ 1 modifier). The result is shown in Figure 03 
In this case, we have integrated out to the virial radius 
of each cluster. We show the variance in 200 stacking 
realizations of the simulated light cone, indicating the 



increase in variation at high flux, as expected due to 
rare, very massive objects projected into some fraction 
of the images. 

We are also interested in the redshift distribution of the 
objects, which can be compared to analytic estimates. 
The result of this analysis is shown in Figure [Sj In each 
case, we have taken the mean value and plotted it as a 
solid line, and the 90% variance in 200 light cone realiza- 
tions as dotted lines. The black lines indicate all clusters 
with total masses above 10 14 M Q in the simulation which 
are in the projected field of the survey. Blue and red 
lines are for clusters above higher mass cutoffs, blue for 
M> 3.0 x 1O 14 M and red for M> 5.0 x 1O 14 M . These 
give a rough indication of the expected redshift distribu- 
tion of identified clusters in upcoming surveys. 

For one projected light cone image, we show the Y 
vs M relationship in Figure [7J In this case, Y is cal- 
culated by integrating the value of Compton y in the 
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Fig. 5. — Plot of number of objects per square degree in the 
simulated survey image as a function of integrated Compton y 
parameter. Dotted lines arc 90% variance as calculated from 200 
stacking realizations of the light cone. 




Redshift 

Fig. 6. — Angular density of clusters in the light cone images as 
a function of redshift. Black solid line is the redshift distribution of 
clusters with M> 10 14 M Q . Blue solid line is for M> 3. Ox 1O 14 M , 
red is for M> 5.0 X 10 14 M Q . Dotted lines are 90% variance of 200 
independent stacking realizations of the light cone. 

image out to each cluster's ■projected virial radius. The 
value of Y is corrected for redshift since it depends on 
E(z)~ 2 / 3 (see iNagal I2006T ) due to the cosmological de- 
pendence of the cluster M-T relation, and the angular 
scale is converted to Mpc through use of the value of 
Da for each cluster. Contrast this plot with Figure 
for a single realization of the light cone image, where 
the true projected Y is plotted against cluster true mass 
from the simulation. The true value for Y is calculated 
from the projected spherically averaged radial profiles 
of each simulated cluster, but includes the gas out only 
to the virial radius in three dimensions. This true rela- 
tion has ex traordinarily tight scatter, as has been shown 
previ ously (|da Silva et all l200fl iMotl et alJ [2005b iNagail 
2006). The difference in the two plots is effectively the 
difference between the cluster's true integrated SZE and 
the SZE integrated in a cylinder with radius equal to 
the cluster's virial radius. For a narrow mass bin around 
3.OxlO 14 M0, the median bias in Y due to projection 
is 79%, and the scatter is +32%/-16% about that me- 
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Fig. 7. — Projected integrated Y value plotted against cluster 
mass. Y is integrated from the simulated light cone survey images 
from the center of each cluster out to the projected virial radius. Y 
is converted from angular units using the angular diameter distance 
appropriate to the redshift of the matched simulated cluster. Y is 
also scaled with E(z) — 2 / 3 to account for the cosmological variation 
of the mass scaling relation. 

dian, a significant increase over the scatter in the true 
Y-M relationship. Note that some clusters in Figure [7] 
appear to have lower value of Y in projection than the 
true value for that cluster. These are the clusters which 
lie near the edge of the simulated survey image, and the 
cluster extends beyond the image edge. As identified by 
other studies, errors in extracting the correct value of Y 
should dominate the error budget fo r this relation (e.g., 
IWhite et all l2002t iMelin et alJ l2006h . These figures il- 
lustrate the difficulty in accurately estima ting Y from 
observ ations. A similar result is noted by IWhite et al.1 
(2002) using data from a smoothed particle hydrodynam- 
ics (SPH) simulation. We have corrected the primary 
object in the source region for cosmological evolution in 
the M-T relation, so the remaining scatter results from 
line of sight projection effects. This means that there 
are secondary bound objects (and unbound gas) which 
are projected into the cluster's source region as will be 
described in the next section. 

4.4. Confusion Problems 

In this section we address the confusion resulting from 
clusters, groups, and smaller mass halos, as opposed to 
confusion resulting from radio and millimeter wave point 
sources in the cluster fields. This work also includes the 
additional flux contributed by unbound gas, which con- 
tributes at some level in the real universe, but has been 
ignored in most simulations of SZE surveys. While there 
are several definitions of confusion in the literature, even 
for SZE surveys, we have chosen to define cluster confu- 
sion as the number of true cluster- or group-mass objects 
in the simulation whose centers are within the source 
region in projection. Since upcoming SZE surveys are 
unlikely to detect cluster gas out to the virial radius, 
we have chosen a smaller radius (7"5oo) as the source re- 
gion. For each of the upcoming surveys, we have defined 
the source region as a radius of a full beam diameter, 
presuming that if two objects were imaged by two sep- 
arate non-overlapping beams, that there is a possibility 
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Fig. 8. — Integrated Compton y vs Total mass relationship ex- 
tracted directly from the simulation data. Y is the projected inte- 
grated SZE y parameter from each cluster out the the virial radius, 
total mass is from the simulation grid for each object. Y is also 
scaled with E(z) -2 / 3 to account for the cosmological variation of 
the mass scaling relation. 

they would not be confused. This also presumes that the 
secondary object is bright enough to be detected on its 
own. One c an also define confu sion as an error in re- 
covered flux (|Holder et al.ll2007t ) from sources found by 
progressive matched filtering at different angular scales 
corresponding to variation of the angular size of clusters 
as a function of redshift (e.g.. iMelin et al.ll2006f) . 

Whether this type of confusion can be mitigated de- 
pends on a variety of factors, including the mass of the 
additional secondary objects in the source region, red- 
shift distribution and angular scale of the objects and 
the observing beam and multiple wavelength identfica- 
tion of the objects (e.g., optical, X-ray, lensing). We 
endeavor here then simply to characterize the amount 
of said confusion, leaving the mitigation of this problem 
to future work. Figure [9] shows for each of the surveys 
considered a histogram of the number of objects above 
10 14 M Q in the source regions (as defined above) of all 
clusters with M > 3.0 x 1O 14 M . These mass limits 
are chosen somewhat arbitrarily, though for the ground- 
based surveys, 3.0 x 10 14 M Q is close to the expected mass 
limit. 10 is chosen for the contaminating objects 

on the expectation that one or more object above that 
mass in the source region should lead to a significant bias 
in the SZE flux from that expected from a single clus- 
ter in the detectable mass range. As in Section 4.3, we 
have identified each cluster by its projected image plane 
position associated with its true three-dimensional loca- 
tion in the simulation after stacking, ft is clear that for 
Planck, there are a high number of objects per source 
region above this relatively high mass. In fact, roughly 
40% of source regions have only one object above this 
mass cut, while nearly 60% have more than one. Con- 
trast this result with that of APEX/SPT, where 95% of 
source regions have only one massive cluster projected 
into them. This result is clearly correlated directly with 
beam size, but also depends weakly on sensitivity, since 
a deeper survey will lead to larger source regions (when 
limited by signal-to-noise) and will increase the chance of 
objects overlapping into the source, and can also increase 
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Fig. 9. — Histogram of number of simulated clusters with M 
> 1O 14 M0 per source region for each of survey modified light cones. 
Source region is defined as inside r'500 projected for the solid line, 
and as within 1 beam diameter (at 144 GHz) distance for each of 
the survey instruments listed. Cluster identifications for all 200 
light cones are tallied and used in the fraction. 

"bridging" between sources. Also, it is likely that SZE 
experiments will sometimes detect a source where there 
is no cluster in the mass range at that location (spurious 
detections). Therefore our estimates shown in Figure [5] 
may be optimistic. 

4.5. Contribution of Unresolved Halos and Unbound 

Gas 

One possibly significant difference between this and 
other similar simulation-based studies is the inclusion 
of adiabatic hydrodynamics in addition to N-body dy- 
namics in the calculation. This results in several advan- 
tages over N-body calculations that include the effects of 
baryons in the post-processing phase. The first is that 
our clusters need not be in hydrostatic equilibrium (HSE) 
which is a standard assumption in dark matter-only sim- 
ulations which have been post-processed. Both simula- 
tions and observations indicate that hydrostatic equilib- 
rium is not a safe assumption for many cl usters (see e.g., 
iRasia et al.ll2005UMarkevitch et alJl2002ft . A significant 
amount of scatter in cluster obser yables results from the 
disequilibrium caused by mergers (iRoettiger et al.lll996l : 
iRicker fc Sarazinl 120011 : [Randall et al.ll2002f) . This scat- 
ter is absent in N-body + HSE type studies, and is nat- 
urally included in our work. Secondly, our simulations 
include baryons which are outside virialized objects, in- 
cluding gas in filaments and voids. The SZE in particular 
is sensitive to this additional gas, since the effect is only 
linear in the gas density and is relatively redshift inde- 
pendent. Thus any gas along the line of sight contributes 
to the SZE integral, and is not diminished by distance. 

While several authors have noted that the angular 
power in the S ZE from the clust er subtracted field is 
small (see, e.g.. lHolder et al.ll2007f ). it is not necessarily 
true that the total SZE flux (or decrement) from unre- 
solved halos and unbound gas is negligible. Figure [TU] 
shows an image of the full field of light cone in projected 
Compton y parameter next to an image of the field where 
clusters with M > 5.0 x 10 13 M Q have been removed from 
the field. We show in Figure [Til the angular power from 
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Fig. 10. — Left panel: 100 deg 2 projected light cone image of the Compton y-parameter from a 512 3 Enzo simulation of a (512/i — 1 
Mpc) 3 volume with 7 dynamic levels of refinement. Light cone includes tiles at 27 discrete redshift intervals between z=3 and z=0.1. Right 
panel: Same image as left panel, but with clusters with M > 5 X 10 13 Mq cut from the data. Roughly one third of the total flux in the 
image comes from the objects that remain after the removal of the massive clusters, including poor groups and filaments. We predict that 
such observations could provide the first detection of the WHIM gas over large sky areas. 

(WHIM) phase. Figure [T2l shows a histogram for 200 in- 
dependent light cone realizations from our simulations of 
the ratio of the total SZE flux (or integrated Compton y 
parameter) in the 100 deg 2 field from only gas within the 
virial radius of those clusters to that of all other gas in the 
field. Thus, we predict that upcoming SZE instruments 
are the only near-term telescopes that will possibly be 
capable of detecting WHI M over large sky areas. This 
result is consistent with iHernandez-Monteagudo et al.l 
(2006), who have performed a similar study with a fixed 
grid N-body/hydro simulation with considerably lower 
peak resolution (195 h _1 kpc) than our work. 

It is important to note that in an adiabatic simu- 
lation, gas fractions are relatively constant with clus- 
ter/group mass. In the real universe, as well as in 
mo re realistic simulation s such as those we performed 
for lHallman et all (|2006f ). the ICM gas radiatively cools 
and forms stars, lowering the gas fraction of the cluster, 
and effectively attenuating its SZE signal. This effect 
is also cluster-mass dependent, lower mass objects have 
comparatively lower gas fractions. Within a simulation 
with radiative cooling and star formation prescriptions, 
we find gas fractions 30-50% lower in clusters with M 
w 1O 14 M than the average value for clusters with M 
> 3 x 10 14 M(D. Additionally, r ecent observational stud- 
ies (e.g.. iMcCarthv et al.ll2007T l point out that gas frac- 
tions deduced from X-ray data decrease with cluster tem- 
perature (therefore with mass), by of order 50% for 1-2 
keV clusters from a flat value above ~4 keV. Therefore 
we expect our result here to be an upper limit on the 
amount of flux from low mass objects and WHIM gas, 
possibly above the true value by as much as a factor of 
two. Our future light cone simulations will be run with 
non-gravitational physics, and we will explore the effect 
on the SZE background. 
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Fig. 11. — Angular power from 200 full light cone images 
(solid line) compared to angular power from images where M 
> 5 X 1O 13 M0 halos are subtracted (dashed line). Dotted lines 
indicate 90% variance range for the 200 independent stacking real- 
izations of the light cone images. 



the SZE in the cluster subtracted images compared to 
the power due to the full SZE image. At small scales 
(£ > 5000), the power in the cluster subtracted image is 
more than an order of magnitude lower than that of the 
full image. At larger scales, however, the difference nar- 
rows, and in fact the 90% variance overlaps for the two 
in some regions. This is likely due partly to incomplete 
cluster subtraction in the image. 

A significant result from this analysis is that roughly 
one third of the SZE flux in the image comes from ob- 
jects with M < 5.0 x 1O 13 M and filamentary structures 
made up of gas in the Warm-Hot Intergalactic Medium 
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Fig. 12. — Histogram of the ratio of total flux in the SZE y 
parameter images from clusters with M > 5 X 1O 13 M0 and from 
images with clusters subtracted. The histogram is generated from 
200 independent realizations of the light cone using the same sim- 
ulation. Roughly two thirds of the flux in the image comes from 
clusters with M> 5 X 1O 13 M0, and the other third comes from the 
WHIM and poor groups. 

4.6. Contribution per Source of Cluster Subtracted 
Images 

Since unresolved halos (M < 5 x 1O 13 M0) and unbound 
gas in this simulation clearly contribute flux to the im- 
age, it is important to ask how much additional flux per 
source is added. This extra flux is a bias, in that it is 
always additive. Therefore it should boost the emission 
of all clusters in the field by some amount which may 
vary from cluster to cluster, adding both bias and scat- 
ter to the cluster SZE observable. This effect is critical 
to understand, since photometric accuracy of the SZE in 
clusters is key to calibrating a Y-M relationship which 
will be useful in determining cosmological paramteters. 
The precision of the calibration of the Y-M relati onship 
depends strongly on its scatter (|Melin et al.ll2006l ). 

Figure [13] shows the ratio of the integrated Compton 
y parameter inside a cylinder of radius r$oo around each 
cluster above 1O 14 M0 from the cluster subtracted image 
to the value for the full image. The plot shows values 
for clusters in a single light cone image, but there is very 
little variation in the mean, median and scatter across 
all the realizations. The large scatter is partly a result of 
incomplete subtraction of the cluster SZE signature from 
the image, but also results from variations in large scale 
structure in the various source regions. What is clear 
from the plots is that there is a systematic bias in inte- 
grated Compton y resulting from low mass objects and 
unbound gas. The mean value of this ratio is 16.3lg4%. 
The la scatter is ±30-40%. It is unclear that this bias 
is reducible, since identifying objects of very low mass 
in source fields, particularly at high redshift, seems pro- 
hibitively difficult. Additionally, making direct observa- 
tion of filamentary gas, particularly to locate its position 
on the sky, has been nearly impossible. 
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Fig. 13. — Ratio of integrated Compton y inside a cylinder 
projected to rsoo from the cluster-subtracted image to integrated 
Compton y from the same projected cylinder in the full light cone 
image. Includes all clusters projected into the survey image which 
have M > 10 14 M Q . 

There will be some variation from survey to survey 
in this additional flux. Depending on how much of the 
cluster's radial extent is sampled, the mean value for all 
clusters will change. The scatter though, is quite large, 
which means that accounting for this flux is not as sim- 
ple as removing a uniform background from each cluster's 
SZE signal. In our study we also performed the identi- 
cal analysis for the clusters, assuming detection out to 
^2500) and found no change in the scatter, though the 
mean value of the additional flux dropped to roughly 8- 
10%. As discussed in previous sections, a more realistic 
modeling of the heating and cooling in the ICM should 
result in a reduction in this additional flux (or decre- 
ment) by as much as 30-50%. Even with a reduction of 
this size, we still expect it to be a few percent to 10% 
effect with of order ±30-40% scatter, creating challenges 
for a percent-level calibration of the Y-M relationship. 

5. DISCUSSION AND SUMMARY 

In this study, we have taken an important step miss- 
ing from previous work in the literature on characterizing 
selection functions of SZE surveys. In earlier work, inves- 
tigators have attempted to evaluate methods of remov- 
ing contamination of the galaxy cluster SZE signal due 
to the CMB and other point sources and backgrounds. 
However, they have included only the baryons present 
in clusters and groups, and artificially inserted the gas 
in hydrostatic equilibrium with the dark matter distri- 
bution from an N-body only simulation. Here, we have 
explored the often neglected contribution of gas in low 
mass halos and unbound filamentary gas in aggregate 
to determine the effect on the cluster signal at a single 
frequency. 

The presence of gas outside the cluster virial radii in 
low density structures such as filaments is potentially an 
important contributor to the cluster SZE signal. This gas 
is completely absent in non-hydrodynamic treatments, 
but appears naturally in our calculation. Since the SZE 
does not diminish with distance and results from a line 
of sight integral of the gas pressure, the sum of all the 
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flux from unbound gas could be a significant contributor 
to the total flux in any cluster detection. Additionally, 
our simulation includes the cluster gas in the full array of 
dynamically active states. Though the integrated SZE is 
less sensitive to cluster-cluster merging than an X-ray ob- 
servation would be, merging contributes non-zero scatter 
to the Y-M relationship. Thus, for all these reasons it is 
critically important to self-consistently include baryons 
in numerical simulations in order to properly simulate 
sky surveys. 

We have shown that on 100 square degree patches of 
the sky at I < 2000 or so the deduced power can be dif- 
ferent by factors of 5-8. This indicates that the power 
spectrum can be quite different from one area of the sky 
to another, and clearly requires much larger angular ar- 
eas to be well constrained. The effect of cosmic variance 
does not become very small until £ > a few thousand. 

We have shown that projection effects can create a 
large bias and additional scatter in the value of Y 
measured for clusters of galaxies. For clusters of M 
«3.OxlO 14 M , the median bias in Y due to projection 
is 79%, and the scatter is +32%/-16% about that me- 
dian, a significant increase over the scatter in the true 
Y-M relationship. Additionally, we have shown that the 
contribution of low-mass unresolved halos and unbound 
gas to the flux (or decrement) of identified sources can 
be significant in some cases, and certainly varies widely 
from source to source. We find that there is a contri- 
bution from gas outside clusters of 16. 3+™% P er object 
on average for upcoming surveys. This indicates both 
a bias and an additional source of scatter in the deter- 
mination of the true SZE signal from any given cluster. 
As identified by other studies, errors in extracting the 
correct value of Y should dominate the error budget for 
the Y-M relation. This effect is critical to understand, 
since photometric accuracy of the SZE in clusters is key 
to calibrating a Y-M relationship which can be useful for 
the precision determination of cosmological parameters. 
While the intrinsic Y-M relation has very small scatter, 
what matters in practice is the ability to determine the 
value of Y accurately. The precision of the calibration of 
the Y-M relationship depends strongly on its scatter. 

We also show results from an analysis of the source con- 
fusion for each instrument based on how many massive 
(M > 10 14 M ) clusters lie within an identified source 
region (within a radius of rsoo or a beam diameter at 144 
GHz). It appears that pure cluster/group confusion in 
these surveys will be a significant problem, particularly 
for Planck Surveyor, but also to a lesser extent for the 
other single-dish surveys. Smaller beam surveys (~1') 
have more than one massive cluster within a beam diam- 
eter 5-10% of the time, and a larger beam experiment like 
Planck has multiple clusters per beam 60% of the time. 
We may have slightly overestimated the problem, since 
the use of higher resolution (shorter wavelength) bands 
on some of the survey instruments will help to alleviate 
this issue. On the other hand, we have not accounted for 
spurious detections which may result in a field with real 
backgrounds and instrument noise. Whether this type 
of confusion can be mitigated depends on a variety of 
factors, including the mass of the additional secondary 
objects in the source region, the redshift distribution and 
angular scale of the objects and the observing beam, and 
multiple wavelength identfication of the objects (e.g., us- 



ing optical, X-ray, and lensing measurements). 

This study uses a large volume, high resolution adi- 
abatic simulation, which serves as a template for more 
complex runs involving additional non-gravitational 
physics that is likely important to accurately modeling 
cluster surveys. There remain important differences be- 
tween simulation outputs and cluster observations, par- 
ticularly for lower mass clusters, which indicates the need 
for a better understanding of the details of baryonic clus- 
ter physics. In addition, deviations from isothermality 
and hydrostatic equilibrium in the cluster gas can have a 
strong impact on both the observable and derived prop- 
erties of clusters. It is also important to note that there is 
some dependence of the SZE signal on the details of the 
ICM physics (heating, cooling, conduction, etc.) which 
is not modeled in this work. In future work, we will 
explore the impact of this additional physics, modeled 
self-consistently within the hydrodynamic framework of 
the simulation code, on a selection of SZE clusters from 
surveys. 

This work will be expanded in future papers by a de- 
tailed treatment of the point source confusion and in- 
strumental and observing limitations. This will include 
adding to our synthetic surveys the contribution of the 
CMB, dusty galaxies, AGN, and atmospheric foreground. 
We are now working on modeling these, in particular 
for APEX-SZ and SPT, and will then experiment with 
matched filtering and the use of multiwavelength cover- 
age provided by SPT for example to mitigate confusion 
and remove the atmospheric and CMB signals. Matched 
filtering is a process by which the images are filtered with 
a kernel matched to the presumed size and shape of the 
expected sources in the field. In the case of the SZE, this 
procedure filters out information on larger scales where 
the primary CMB anisotropies will be a source of con- 
fusion. It also maximizes the contrast of the image for 
the objects at that scale. Since the angular scale sub- 
tended by massive clusters is a function of redshift and a 
weak function of cluster mass, spatial filtering will need 
to be done at a varie ty of angular scales to get a complete 
cluster sample as in iMelin et al.l (I2006T). A simil a r type 
of analysis has been performed by ISehgal et al.l (2006) 
for ACT's survey. 

We also are currently performing additional synthetic 
light cone sky surveys at X-ray wavelengths, to take a 
first look at the limitations of current (e.g. XMM-LSS) 
and upcoming (e.g. eRosita) X-ray surveys in extracting 
cosmological parameters. 
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